Isotropic Cosmological Singularities II. 
The Einstein- Vlasov system. 



K.Anguige 
K.P.Tod 

February 7, 2008 



Abstract 

We consider the conformal Einstein equations for massless colli- 
sionless gas cosmologies which admit an isotropic singularity. After 
developing the general theory, we restrict to spatially-homogeneous 
cosmologies. We show that the Cauchy problem for these equations 
is well-posed with data consisting of the limiting particle distribution 
function at the singularity. 



1 Introduction 

In the accompanying paper, (Anguige and Tod 1998, hereafter ATI), we 
considered the Cauchy problem for the conformal Einstein equations with a 
perfect fluid source and with data given at an isotropic singularity. Among 
other results, we were able to show that a perfect fluid cosmology, with a 
polytropic equation of state, for which the Weyl tensor vanished at an initial, 
isotropic singularity was conformally-flat: if the Weyl tensor is zero initially 
then it is always zero. This is a strong result, and one naturally wonders how 
robust it is. In particular, does it continue to hold with other matter models? 



In this, the second paper of the series, we consider the same Cauchy 
problem but with a source given by a different matter model, a massless 
collisionless gas. This is the Einstein- Vlasov system. Our long-term aim is 
to prove existence and uniqueness for this system with data at an isotropic 
singularity. However, in the present paper we only partially achieve this 
aim: we are able to prove what we want only with the restriction of spatial 
homogeneity. This is sufficient to investigate the Weyl tensor problem de- 
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scribed above, but of course the long-term aim remains. 

In the Einstein- Vlasov system (see e.g. Ehlers 1971, Rendall 1992, 1997) 
the matter is determined by the distribution function f(x a ,pb), a positive 
function on the cotangent bundle of space-time, which is interpreted as the 
density of particles at the space-time point labelled by x a which have 4- 
momentum p^. The collisionless assumption is translated as the assumption 
that the distribution function is constant along the geodesic flow, giving 
an equation known in this context as the Vlasov (or Liouville) equation. 
The Einstein equations are written down with a stress-tensor defined as an 
invariant second moment of / integrated over the fibre of the cotangent 
bundle. Thus / is constrained by the condition that this moment exist. We 
will consider only massless matter, which means that / is supported on the 
null-cone at each point. For analytical convenience, we will assume that / 
is supported compactly on the cone and away from the vertex of the cone. 

The massless assumption will imply that the theory has a good trans- 
formation under conformal rescaling of the metric. For example, the dis- 
tribution function is constant along null geodesies, which are conformally- 
invariant. Also the stress-tensor is trace-free, which is the case for radiation 
perfect-fluid where it is known from ATI that the conformal factor can be 
chosen to be smooth at an isotropic singularity. 

We refer to the accompanying paper for the background to conformal 
rescaling. Recall only that the physical metric g a b of space-time M is related 
by rescaling to an unphysical metric g a b on a larger space-time M according 
to: 

9ab = Q 2 9ab 

where the conformal factor f2 vanishes at a smooth surface £ in the extended 
space-time M. As remarked above, the massless assumption will enable us 
to suppose that f2 is smooth at S. 

After establishing the conformal transformation properties of the mat- 
ter model, we shall seek Bondi-type expansions of space-time quantities as 
power-series in the conformal factor, which becomes a good time-coordinate 
in the unphysical (rescaled) space-time. In this way, we are able to iso- 
late the degrees of freedom of the problem, or equivalently to identify the 
Cauchy data. The data are very different from the perfect-fluid case. We 
need the first and second fundamental forms of the singularity surface S 
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and the initial distribution function / . These turn out to have very strong 
constraints among them, so strong in fact that the initial distribution func- 
tion f° uniquely determines the first and second fundamental forms of the 
singularity surface. The initial distribution function is subject to a single, 
integral constraint which can be interpreted as a 'zero drift' condition. 

Recall from ATI that in the perfect-fluid case, the second fundamental 
form of S is necessarily zero and the first fundamental form is free data. 
There is no free data for the matter (apart from the equation of state). 
One can say aphoristically that the metric determines the matter. From 
the previous paragraph, we see that circumstances are quite different in the 
Einstein- Vlasov case where the matter determines the metric. 

Because the second fundamental form of the singularity necessarily van- 
ishes for perfect fluid, the condition of initially-vanishing Weyl tensor forces 
the physical cosmology to be FRW, so that, as we noted above, the Weyl 
tensor is always zero if it is initially zero (ATI). With Einstein- Vlasov the 
situation is rather different. The second fundamental form is determined by 
/° and its vanishing is an integral constraint on /°, related to the vanishing 
of a third moment. It is quite possible, as we shall see in section 6.6, that 
the lower moments of f° are such that the initial Weyl tensor is zero but 
there are higher moments which allow the later evolution to diverge from 
FRW, so that non-zero Weyl tensor emerges. 

Once we have identified the data, we turn to the problem of existence 
and uniqueness of solutions. For this, we restrict to spatially- homogeneous 
cosmologies. Then the Einstein equations become a system of ordinary 
differential equations with a singularity in the 'time' and a source involv- 
ing integrals over the distribution function /, while the Liouville equation 
remains a partial differential equation for /. We prove existence and unique- 
ness for this coupled system by a method based on results of Rendall (1994) 
and Rendall and Schmidt(1991). 

The contents of the paper are as follows. We begin in section 2 with 
a review of relativistic kinetic theory and then in section 3 consider the 
problem of conformally-rescaling the massless Einstein- Vlasov (henceforth 
EV) system. The Bondi expansions are begun in section 4, which includes 
the proof that the constraint relating f° and the metric a^- of S determines 
aij uniquely. In section 5, we fix the conformal gauge for the EV system 
which enables us to simplify the expressions found for the data in section 4 
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by eliminating some gauge quantities. Section 6 is the heart of the paper. 
We adapt the expansions of section 4 to a spatially-homogeneous cosmol- 
ogy, identify the data and then prove, in Theorem 6.1, that solutions exist 
and are unique given suitable data. We consider what data lead to FRW 
cosmologies and show in particular that there are Bianchi-type V solutions 
in which the Weyl tensor is zero initially but not for all time. 



2 Relativistic kinetic theory 

In kinetic theory, the matter content of spacetime M is taken to consist of 
a collection of particles which move on geodesies between collisions. The 
motion of a freely falling particle is determined by the Lagrangian: 

L(x\x d )= X -g ah ± a x h (1) 



since the Euler-Lagrange equations are equivalent to the (affine) geodesic 

8L 

dx a 



equation. Writing p a = g a bX b = Mh = canonical momentum, we get the 



Hamiltonian 

H(x c , Pd ) = l -g ab Pa p b (2) 

and thus we obtain canonical equations for the geodesies. H is constant 
along each particle world-line, p a is the 4-momentum, and m = \/p a p a is the 
mass of the test particle. To take into account all possible particle positions 
and momenta, we must consider the one-particle phase space P . This is 
the collection of instantaneous states (x a ,g a bx b ) of a single particle and is 
the future causal part of the spacetime cotangent bundle. 

States having a given mass m > form the seven-dimensional subman- 
ifold P m of P. On P m we take (x a ,pi) as local coordinates, po being 
determined by the equation g ab p a Pb = m2 , an d the requirement that p a be 
future directed (conventionally, a,b,c, .. = — 3 while k, .. = 1 — 3). The 
free-fall trajectories define on P a congruence of curves, along which H is 
constant. The geodesic spray £ is a vector field along these curves defined, 
in local coordinates, by 

ab 8 1 dg ab d 
C = 9 Pa d^-2 PaPb ^dp- c (3) 
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The cotangent space to M at x is a flat Lorentzian manifold, and on the 
submanifold P m (x) there exists an invariant volume measure Tr m given by 

vr m = ^f^V (4) 

The distribution of particles and momenta is described by a scalar func- 
tion / = f(x a ,pt)) on P, and the central result of relativistic kinetic theory 
is the equation 

£(/) = 9 (5) 

where g describes the density of particle collisions. Hence the condition 
that / represents a collisionless gas is simply 

£(/) = (6) 

and this is known as the Vlasov equation. Note that / satisfies the Vlasov 
equation if and only if it is constant along geodesies. 

The stress-energy-momentum tensor due to particles of mass m is given 

by 



1 ab 



(x) = / fPaPb^m (7) 

JPm(x) 



'Pm{x) 

and if the Vlasov equation is satisfied then 

V a T™ b = (8) 

If the gravitational field is due solely to particles of a single mass m, rep- 
resented by /, then the coupled Einstein- Vlasov equations, for the met- 
ric g a b and the particle distribution / are 

G ab = 8TT fPaPb^m (9) 

J Pm 

W) = (10) 



3 Conformally rescaling the massless EV equations 

Now we shall restrict to massless, collisionless particles and investigate the 
transformation properties of the Einstein- Vlasov system under conformal 
rescaling. For a future-pointing, null co- vector p a , the component p$, either 
coordinate component or orthonormal-frame component, is determined by 
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the spatial components pi. Thus we write f{x a ,pi) for the distribution 
function for massless particles on Pq in a spacetime (M,g a [>). Suppose that 
g a b and / satisfy tilded versions of (9)-(10). Let g ao be the metric defined 
by 

9ab = ft 2 9ab (11) 

Under this transformation, null geodesies of g a \, become null geodesies of g a b- 
Thus, since / is constant along the geodesic flow in (M,g a i,), f = f will be 
constant along the geodesic flow in (M,g ao ). In fact one calculates that 

c~J = n- 2 c g f (12) 

and it follows that / satisfies the Vlasov equation in M. 
If we now define 

{- 9 )z 112 



Tab = J fPaPb d P ( 13 ) 

then the conservation equation X7 a T ao = holds in M, just as the tilded ver- 
sion holds in the physical spacetime. This follows either from the unphysical 
Vlasov equation, or the relation 

f ab = ^T ab (14) 

together with = 0, which in turn follows from m = 0. 

With the conformal transformation of the Ricci tensor taken e.g. from 
ATI, we can now write down the conformal EV equations for g^ and / as 

R ah = 2v a v fe log n - 2v a log nv b log n 

+ <7 a6 (Dlogfi + 2V c log^V c log^) + ^y fPaPb { ~ 9) p0 d 3 p (15) 



and 

Note that (15) implies 
since T n a = 0. 



= o (16) 



nn = -im (17) 

6 v ' 



Now, as Vl is smooth in M and satisfies the regular, linear wave equation 
(17) in M, one must have that O has a smooth extension onto M (see Racke 
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1992). It follows that Q is a smooth function in M. 



In the sequel it will be important that be a good coordinate in M, 
and thus for we make the following assumption: 

Assumption 3.1 The conformal factor $7 is such that V a ^ ^ at S. 

Thus corresponds to the time function Z of ATI in the case of a radi- 
ation fluid, and for convenience we will henceforth write £1 = Z. 



4 The initial data: Bondi expansions near S. 

The singular nature of the field equation (15) imposes constraints on the 
data at the singularity £ . To see what these constraints are we shall seek 
'Bondi' expansions of the field variables in powers of Z near £ (analogous to 
the corresponding expansions near future- null-infinity, whence the name). 
That is, we write 

hij = h% + Zh\j + ... 

etc, and compare coefficients in (15)-(16). To aid calculation we use coordi- 
nates x a , with x° = Z and the x i , (i = 1, 2, 3) taken to be comoving along 
the integral curves of V a Z. In these coordinates the line element in M can 
be written 

1 

V 2 

with h^ positive definite and V 2 = g ab Z^ a Z^. 



ds 2 = —dZ 2 - h ij dx i dx j (18) 



The Vlasov equation for / = f(x a ,pi) takes the form 

V2P °% ~ kijpi §xl " \^ V2 ^ 2 ~ (^ hmn )P m Pn}^ = (19) 

where V 2 (po) 2 = h^piPj. 

The Einstein equations are just 

Z 2 R ab - 2ZV a V b Z + 4V a ZV b Z = g ab {ZUZ + V 2 ) + 8nT ab (20) 

4.1 O(l) terms 

Considering the lowest order terms in (20) gives 

87r(Too)° = 3, 8vr(r 0l ) = 0, ^{T %J f = {V°) 2 h%. (21) 
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The second of these is just an integral constraint on /, the distribution 
function at E: 



where h° = det(/i°-) and the first is implied by the third since is trace- 
free. 

It is now natural to ask two questions: 

1. Given a positive, suitably integrable function /°, does there exist a 
3- metric satisfying equation (23)? 

2. Given a positive-definite 3-metric h®j does there exist a positive func- 
tion f° satisfying equation (23)? 

These questions are answered affirmatively by Theorems 4.1 and 4.2 below. 

Theorem 4.1 Let f{x->,Pi) be a smooth, positive function on (7 x R 3 , 
where U is an open subset of M 3 . Suppose that for each x € U 

1. / is compactly supported in p; 

2. / is supported outside some open ball containing p = 0; 

3. / is not identically zero in p. 

Then, given a smooth strictly positive function V{x) on U, there exists a 
unique positive definite 3-metric hij(x) on U satisfying (23). Moreover, the 
metric hij(x) is smooth. 

Proof. The metric is found by a minimisation property. First define the 
function F by 




(22) 



the third gives an implicit relation between f° and h -: 




(23) 



F(6 n ,6 22 ,6 33 ,6 12 ,6 




where 




(25) 
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and B is positive definite. 



Clearly F(B) = F(XB) for A £ R, so we can think of F as a function de- 
fined on a region S 1 in R 6 consisting of points P = (6 11 , b 22 , 6 33 , b 12 , 6 13 , 6 23 ) for 
which B is positive definite and tr£? = 1. 

Now for a couple of Lemmata: 

Lemma 4.1. The function F has a critical point: 

|^ = * = M,3 j>i (26) 
if and only if the following equation holds 

^jtfBtiWftv-ijjgSfa*, (27) 

where (a^) = (6 y ) _1 . 
Proof. Note that 

S(detB) 26ij 



where 6^ = (det-B)ajj, and there is no summation on the rhs. 



Hence 

dF 



iCdetBj-Vfi-L-^. J(p T Bp) l ' 2 f d 3 p 



and the lemma follows. 



Lemma 4.2 If M = (m !J ) is a 3 x 3 matrix with trM = 0, then 



(28) 



+ (detS)- 1 / 6 ^ j (p T B P r 1 / 2 p t p 1 f d 3 p (29) 



B 2 F 

S 2 F = ". n ,., m ij m kl 2( s ><'- 1 W s v-V > (30) 
dWdb™ 



at any critical point of F. 
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Proof. By definition one has 

bij = —6ipq€j rs b P W (31) 

which implies 

Q\) 1 

~ n(Ski) ( €i P ke i rl e ipi e jrk)ti p (32) 



and thus 
d 2 F 



" 2(1 ^ | 6 



-if(-B) ("-((detB)- 2 )2< 1 -*»'i M i a + (detfl)" 1 ^(e^ejrt + e ip ,£j> t )(W r )j 

(33) 

At a critical point we may, without loss of generality, take B = 1% , the 
3x3 identity matrix, for if 6 y solves (27) then 1% solves the same equation 
with / replaced by /, where f(x,p) = f(x, hp) for some linear transforma- 
tion L, and we can work with /. 

Now drop the hat and evaluate (33) at a critical point to give 

d2F - 2 (l-^) 2 (l-^0/ _ I f fPiPiPkPl d 3 J_ g f fPiPj d 3 
d¥W~ \ 4 J (pTp)3/2 aP 12° kl J (p T p) 1 / 2 P 

- ^F(h) ( - hAj + ^{eipktjrl + e ip ie jrk )5 pr ^j j (34) 

so that 

db l Wb kl 4 7 (p T p) 3 / 2 12 7 (p T p) 1 / 2 

~ ^F(h)(-(trM) 2 + 2trM) (35) 
Note that for any 3x3 matrix M there is the following relation: 

trM = ^{(trM) 2 - trM 2 } (36) 
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where M is the matrix of cofactors. Hence (35) becomes 

Now let the eigenvalues of M be Ai,A2,A3 with 
Ai > A2 > 0, and A3 = — (Ai + A2) (replacing M by — M if necessary). 
It follows that 

\p T Mp\ 2 < |A 3 | 2 b T p| 2 (38) 

and equality holds for only one value of p (up to scale) . 
Hence 

5 2 F> -Lf(/ 3 )(Ai-A 2 ) 2 >0 (39) 
Note that we have actually shown that F, defined by 

F(B) = (detB)- 1 / 6 J {p T Bp) 1 / 2 / d 3 p 

is such that 5 2 F > at B = I3. But we have 

F(B) = |detLr 4 / 3 F((L T )- 1 5L- 1 ) 
It follows that 5 2 F > at a critical point B. Hence the lemma. 

Suppose next that (27) holds and put hij = Aa^- where 

Then solves (23) and so there is a one-one relation between critical 
points of F and solutions of (23). 

Returning to the set S on which F is defined, we see easily that it is both 
convex and open in O = M 6 n {B : trB = 1}. Since B is +ve definite we get 
from (36) that (trB) 2 > trB 2 and consequently max{|6 12 |, |6 13 |, |6 23 |} < |, 
so that S is bounded. 

Since S is convex it must just consist of radii from some point in S, in the 
directions spanned by the b 12 , b 13 , b 23 axes, and in the plane given by 

b ll + b 22 + b 33 = 1 

The convexity also implies that S has continuous radius, and thus is home- 
omorphic to an open ball in R 5 . In particular S is contractible: it has the 



11 



homotopy type of a point. By Lemma 3.2 S 2 F > at a critical point, when 
variations are restricted to directions in S, so that any critical point of F is 
a minimum. 

Now F — ► oo on the boundary dS of S (thinking of S as a subset of 0), 
so that j, is continuous on S, with j, = on dS. Now choose a > and 
consider the following subset Q of S: 



Since 5 is bounded and y is continuous, one must have that Q is a closed, 
bounded subset of O, disjoint from dS. Writing G = —F and b = —a now 
gives that oo) is compact. Since the critical points of G are all max- 

ima, it follows from theorem 3.5 of (Milnor 1963) that S has the homotopy 
type of a discrete set of points, one for each critical point of G. But we 
already know that S has the homotopy type of a single point and thus F has 
just one critical point, as required. 

It remains to show that the metric determined by (23) is smooth. To do 
this first write 

dF 

Ha = (40) 

and suppose that B° is the critical point of F at x°, so that Hij(B°, x°) = 0. 
Now, by hypothesis, Hij is smooth in B and x in a neighbourhood of (B°, x°), 
and we have 

™» = (41) 
dh ki db kl dW K ' 

But from Lemma 4.2 we have that 5 2 F > at (B°,x°), and thus (dsH) is 
invertible at this point. From the inverse function theorem it now follows 
that B{x) defined by Hy(B,x) = depends smoothly on x. 

Theorem 4.2 Let hij(x) be a smooth, positive-definite metric on an open 
set U C M 3 , with V a smooth, strictly positive function on U. Then there 
exists a smooth positive function / on U x M 3 satisfying (23), and in fact 
there are many such / for each hij. 

Proof. First note that equation (23) can be written in matrix form as 

(VvasH/fcrW) = / ^ Z H %„ " 3 p (<*) 

where H = (hij). 

Now let 4> : R — > R be a smooth function and n, r2 positive constants, 
chosen so that 
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1. <f>{x) = for \x\ < r\ 

2. <p(x) = for \x\ > r2 

3. <f>{x) > for n < |ar| < r2 
4. 

y <f>(p T p) T 3 

h = iWvW w p 

If we let f(x,p) = V 2 (x)(j)(p T Hp), then / is smooth and solves (42). Clearly 
there is considerable freedom in the choice of <p. 



4.2 O(Z) terms 

Let K a b be the second fundamental form in M, so that Kij = —^Vd z hij. 
A calculation of the O(Z) terms in the space-space components of (20) gives 

(28\8\ - X l \) (K°) kl = D m ( X ^ m ) + \(K°)(h yi (43) 

where K° = (h^K^ 

= 4?r f f°PiPjPkPl ,3 (AA ^ 

^ ~ (^0)2^ J m mn pmPn) S/2 P ^ 

( yo )2 ^oy ((^)-> m p n ) dP (45j 

and -Dj is the covariant derivative operator associated with (h°)ij. 
Clearly the eigenvalues of Xijkl regarded as a linear transformation on 3 x 3 
symmetric matrices are positive. Also, from (23), the trace of this linear 
transformation is 3/2, so that each eigen- value is no greater than 3/2. Thus 
the 'matrix' acting on (K°)ij in (43) is invertible, so that the trace-free part 
of (K°)ij is determined by f° through this equation. The trace is not deter- 
mined but we shall see in section 5 that this can be removed by fixing the 
conformal gauge. 

The O(Z) terms in the space-time and time-time components of (20) 
impose no more constraints on the data at Z = 0. 
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5 Conformal gauge fixing in EV 



Inherent in the definition of an isotropic singularity is the freedom to choose 
a new conformal factor Q and conformal metric g ab by rescaling with a reg- 
ular non-zero function 6 according to fi — > Cl = @Q. We now use this 
freedom to make a conformal gauge choice which should be useful for solv- 
ing the conformal EV equations in full generality 

First recall the definition of V as V 2 = g ab V a ZV b Z, and suppose one has 

g ab = Z 2 g ab = Z 2 g ab (46) 
Write h = Z/Z. Then at £ one has 

V 2 = h 4 V 2 (47) 
Since V a Z ^ at S, it follows that h can be chosen so that V = 1 there. 

Now write g ab = h 2 g ab , for some new function h, with h = 1 at S. If ^ is 
the trace of the second fundamental form of g ab , with K defined similarly, 
then at S one gets 

K =^(K + 3N d V d logh) (48) 

where iV a is the unit normal to S with respect to (? a b. We can get K = 
at S by choosing /i so that 

iV d V d fc = (49) 

there. 



Consider now the following wave equation for h: 

Uh = ^Rh (50) 

with initial data h = 1, N d V d h = -\K. 

It is standard (Racke 1992) that this equation has a unique smooth 
solution h with the given data. If we let Z = Z/h, then g ab = Z 2 g ab , and 
from (17) there follows 

DZ = (51) 
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We can now prove the following: 

Lemma 5.1 If (M,g a b) is a solution of the massless Einstein- Vlasov equa- 
tions, with an isotropic singularity, and Assumption 3.1 holds, then the 
conformal factor Z may be chosen so that V(0) = 1, K(0) = 0, and Z is 
a harmonic function in M. Moreover, these three choices fix the conformal 
factor completely. 

Proof. Suppose = Z 2 g a \ ) = Z 2 g a i, , with Z, Z chosen as above so that 
V(Q) = V(0) = 1 , K(0) = K(0) =0, UZ = DZ = 
Define h by h = Z/Z, so that h(0) = 1 and (N d V d h)(0) = 0. One has 

R = h~ 2 (R- -(Dh)) (52) 

But from (17) R = R = 0, and thus 

Dh = (53) 

Hence h = 1, and Z = Z. 



The initial data set for the conformal EV equations at Z = 0, in the 
gauge of Lemma 5.1, can now be summarised as follows: 

• The initial distribution function /° is free data, subject only to the 
integral constraint (22). 

• The initial 3-metric h®j is determined by f° via (23) 
with V° = 1. 

• The initial second fundamental form Kfj is determined by f° via (43) 
with V° = 1, trK° = 

6 The Cauchy problem 

We have seen that the Einstein- Vlasov system has good behaviour under 
conformal rescaling and that it is possible to specify Cauchy data for the 
rescaled equations at an isotropic singularity. In this section we show that 
if attention is restricted to spacetimes with Bianchi type spatial symmetry, 
then this Cauchy data determines a unique solution of the field equations. 

In this section only we will take the spacetime metric g a t to be of signa- 
ture (- + ++). 
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6.1 Bianchi type spacetimes 



The Bianchi type spacetimes are a class of spatially homogeneous cosmolog- 
ical models having a 3-parameter Lie group G of isometries, transitive on 
space-like 3-surfaces, the surfaces of homogeneity. They are classified into 
nine types according to the nature of the Lie algebra associated with G (see 
Wald 1984 for details). In such spacetimes there exists a cosmic time func- 
tion t and one-form fields (e l ) a such that the metric can be written as 



The (e*) a are left-invariant one-forms, preserved under the action of G, and 
the Bianchi time t is proper-time on the congruence orthogonal to the sur- 
faces of homogeneity. 

When field equations are imposed on a spatially-homogeneous space- 
time, they typically reduce to a system of ordinary differential equations for 
the spatial metric (iij(t), together with some matter equations. 

6.2 Conformal gauge fixing for Bianchi types 

Clearly we may assume that the conformal factor Z for Bianchi types is a 
function only of the cosmic time t. There is a natural choice of this function, 
different from that made in Lemma 5.1, given by the following lemma. 

Lemma 6.1 If it is assumed that the conformal factor satisfies V a Z ^ 
at t = 0, then we may take Z = i 1 / 2 , where t is the Bianchi time. 

Proof. Suppose g a b = Z 2 g a t, with g ab smooth and g a b a Bianchi type 
solution of the EV system. Let V 2 = -g ab V a ZV b Z, with V = 1 at t = 
without loss of generality. Define the positive function g by 



Now let g a b = g 2 g a b , Z = gZ. One then gets V = 1/2, Z = i 1 / 2 , and it 
remains to show that g is smooth. 

One must have y _1 (s) = 1 + h(s) for some smooth h with h(0) = 0. 
Hence h(s) = sr(s) for some smooth r, and 



9ab = -VatVfet + a ij {t){e l ) a {e :) ) b 



(54) 




(55) 
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It follows that g is continuous. Now r(s) = r(0) + sp(s) with p smooth, and 
hence 



It now follows that g is C 1 , and continuing in this way, one shows that g is 
smooth. 

6.3 Evolution equations in M 

We now wish to write the conformal Einstein- Vlasov equations as an evo- 
lution system for the unphysical metric g ab and the particle distribution 
function / on a manifold M = [0, T]xG. We assume that / = f(Z, Vi) , 
where the V{ are components of 4-momentum in the invariant frame (e l ) a . 
Note that the metric g a b in M takes the form 

9ab = -4S7 a ZV b Z + a tJ (Z)(e l Ue j ) b (56) 

with ciij positive definite. 

The Vlasov equation can be written 

^ = dzf = 2(b rs v r v s )- 1 / 2 b ln v k v l C^ (57) 

where = (a^) -1 and Cj fc are the structure constants of the Lie group 

G. The spatial components of the Einstein equations, in the basis {e r ) a are 

dz(aij) = hj (58) 

d z {b ij ) = -b m V m k mn (59) 
dz(kij) = ±C k ck (C r t] a tr + a jr Cl)b ct + AC c kt (C% + a cm b kt C%) 

+2C sk C tc a im a ir b b + -^-—^ J d v 

- |% ~ {^Oij + hij^j (b rs k rs ) + b lq k u k qj (60) 

We wish to remove the second-order pole in Z in (60), so we introduce a 
new independent variable Z; L j according to 

1 J 327T f fViVj 3 



z -' = z{(i^Jw^ dv - a ''i <61) 
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Then (60) becomes 

dzikij) = ±CUC' t] a„ + a 3r Cl t )b kt + 4Cfc(Cj + a cm b kt C%) 

+ ^C^Cl c aj m ai r b kt b cs + — {2Zij — Ikij — aij(b rs k rs )) — —kij(b rs k rs ) + b lq knkqj 

(62) 

and using (57) we get an evolution equation for Z^ : 

Note that if (57)-(59), (62)-(63) are solved, and we define Z^ by the rhs of 
(61) then one gets 



dzi^ij — Zij) — — — (Zij — Zij) 



and hence Zij = Zij , and the definition (61) is recovered. Also, from 
(58)-(59) there follows 

dzioijV" - 5?) = -(a tJ V s - 5i)k sn b kn 

so that, by a Gronwall estimate, if QijbP k = Sf at Z = , then a.ijb> k = 5 k for 
all Z. 

6.4 Constraints and initial data 

The Hamiltonian constraint is 

<s£^ / « 6 ""-» 1/2 ^ = I + (3)R 22 - (is) to*** 

+ (f ) W^ojj + ^(l."fc„) 2 (64) 

where ( 3 )i? is the Ricci scalar of a^. 
The momentum constraint is 

/ d " v = Z2 ^% C U + b sl hiC r sr ) (65) 
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Suppose that (57)- (59), (62)-(63) are satisfied, and write (64)-(65) as C = 
0, Q = respectively. Then one calculates that 

d z (Z\det a)C) = , d z Ci = ~(b rs k rs )Ci 

Hence C = and, by a Gronwall estimate, Q = if Q = initially. 
To satisfy the momentum constraint at Z = one must have 

y /V d 3 v = (66) 

which is the counterpart of (22), and if the evolution equations are to be 
satisfied, then from (61) 

n _ 32^ f f°v iVj « 



a H = 71 n\77? / 7 run m /? d v ( 6T ) 

J (det a°) 1 /2 _/ (^0^1/2 

which is the counterpart of (23) (recall here V° = 1/2). If (67) is satisfied, 
then by taking the trace we see that the Hamiltonian constraint is also sat- 
isfied initially. 

To determine k^, Zf,, note that from (62)-(63) 



and 



Z%-k%-\a%({b°yX s ) = V (68) 

+ (det a°)V2 ( ( j ( j mn + [ ' { ' rs) J {v T Wvfn d V ~ U 

(69) 

Now eliminate to get 

2k° + tr&%° WA^V- f | ^ C k (b°) ln I 9f ° v jW!l d \ _ n 
+ trfc (a„j (fc j X «e/+ (det o0) i/2 C n P (ft J y ^ ^T^) d W " U 

(70) 

where 
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Taking the trace of (70) gives 

-^ 0+ (d^ c >°»'"/^' rf3 " (72) 

and then (66) implies that tr/c° = 0. Equation (70) now becomes 

(Xijef-2aiay(kr f = - (^1/2 ^°)^ / jj^v^vt d*v (73) 

which corresponds to (43). By the eigenvalue properties of XijH noted in 
section 4.3, k^- is uniquely determined by /°, and Zfj = k^ by (68). 

In summary, the free data at Z = consists of just f°(vi) satisfying the 
integral constraint (66), with a?-, kfj = Zfj being determined by (67), (73) 
respectively. Given such data the constraints are preserved by the evolution, 
and hence the task is just to solve (57)-(59), (62)-(63). 

6.5 Solving the evolution equations 

We will now use Theorem 1 of (Rendall and Schmidt 1991), together with an 
iteration technique due to (Rendall 1994), to solve the conformal Einstein- 
Vlasov equations near Z = 0. The result, after a series of lemmas, is Theo- 
rem 6.1 below. 

Choose /° G C^°(E 3 ) , /° > 0, with f° supported outside some open 
ball centred on the origin, and let a°-, k^ = Zfj be determined as above. 
Define 

rriij = kij — k^j , Wij = Zij — Z^ 
so that m% = w%, = 0. 

The evolution equations can now be written as follows 

d z f = 2(b rs v r v s )-^ 2 b ln v k v l C^ (74) 

dz{aij) = k^ + mij (75) 
d z (b^) = -b in V m (k mn + m mn ) (76) 
dzirmj) = 4C c kc (C r tj a ir + a jr C^)b kt + 4C c kl (C^ + a cm b kt C%) 



20 



+ P t , mn (k° mn + m mn ) + ±{- a %(b°r n m mn + 2 Wij - 2m t] ] (77) 
dz(wij) = + S$ + l| - Wij - rriij 

i I 1 ^ 7r /" f° v i v j v eVf J3„,\ffifl\enfip\fm m flfl\ef /iO\mn m \l 

+ l v (deta0)i/ 2 y (^^3/2 du J « 6 ) (&)'(&) m mn )j 

(78) 

where 

P ™(*) = ( l { aij V m ¥ m (kl q + m w ) - 6 m "(fe0, + mji )}(st) (79) 
Jo 

and 

ijni.rui.sv / 1,0 i „ „ \ [ df ViVjV k ViV r V s 3 



^ W - l*mn + "W) ^ (a)1/2 |2C pr 6 y ^ {vTbv)2 d V 

+(k° rs + m rs ) ^ r b gs (b en bf m -b e fb mn ) J fVi ^f/f Q d " v 
+ (/^^ d3v )( 6rS(6Cn6/n - 6C/6mn) 



— b er b ns b^ m — b en b^ r b ms + b mn b er b^ s + b e ^b mr b ns ^j 
where a =det(ay). 



\st) ds (81) 



Suppose now that f° is supported outside a ball of radius B\ centred 
on the origin. Let r = {v T v) 1 / 2 , and let 4>{vi) be a smooth function such 
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that (j) = 1 on r > B 2 , B 2 < Bi, <j> = on r < B 3 , B 3 < B 2 , and < <j) < 
1 elsewhere. Replace the Vlasov equation (57) by 

|£ = 2 ( f>(b rs v r v s )-^% ln v k v l C^ (82) 

Equations (77)- (78) can be written 

d z (u) + ^Nu = G(a,b,u,f) (83) 

where u stands for m, w, and N is a constant matrix containing just o°-, (b )^ . 
A" has desirable properties with regard to application of the theorem of Ken- 
dall and Schmidt (1991), given by the following lemma: 

Lemma 6.2 All the eigenvalues of A" have positive real part and A" is 
diagonalisable. 

Proof. The eigenvalue equation for A" is 

2rriij + maf-j — 2wij = Xrriij (84) 

rriij + -malj + ^ - Xijkim kl = Xwij (85) 

where rriij, u;^ are components of an eigenvector, m = (b ) 1 ! rriij, and A is 
the eigen-value. 

Hence Wij = ^(2 — A)mjj + |ma°-, and substituting in (85) gives 

(A 2 - 3A + 4)m„ - (A - 2)moJ- = 2 XijH m fc! (86) 

We know from (67) that (b°)^Xijki = 5^1, and so taking the trace of (86) 
gives 

(A - 3) 2 m = (87) 

If m ^ then A = 3 and an eigenvector is = rriij = \ma^ . If A / 3 
then m = and (86) becomes 

(A 2 - 3A + 4)m ij = 2m kl X ij kl (88) 

and we must consider eigenvalues of X l \i, which we may regard as a 9 x 
9 symmetric matrix, with 9 real eigenvalues and 9 linearly independent eigen- 
vectors. So suppose 

X ij kl m kl = fim ij (89) 
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Then rriij = rriji and we may work in an invariant basis where a ^ 
Sij, (rriij) = diag(mi, 777.2, m 3 ). Clearly fj, must be positive since 

Xijkim ij m kl > 

It follows that 



n\rrii\ = 167T 



fO 

vf(vfmi + v\rti2 + v^ms) d?v 



{v T vf/ 2 



[ f°v 2 

< 16-7T / - 7^ (vi\mi\ + v\\m 2 \ + u 3 |m 3 |) d 3 v 
J [v 1 v) A i z 

and hence 

\ mi \ - 167F / / T \l/2 ^ll mi l + V 2\ m 2\ +^ 3 l m 3|) d 3 V = -^2\ mi \ 
i=l ^ ^ V V > i=l 

by (67), so that fi < ^. 

Going back to (88) one gets that A = ±(3 ± [8/x - 7] 1/2 ) with < \i < 
\. Thus A has positive real part, and there exist two X's for each of the 
nine linearly independent rriij ■ One has Wij = ^(2 — A)my, and this gives 
18 linearly independent eigenvectors of N. 

Lemma 6.2 gives that there exists an 18x18 matrix L such that equations 
(83) may be written 

dz{y a ) + ^Ky a = H a (a, b, y, f) (90) 

where y = Lu, H(a, b, y, f) = LG(a, b, L~ l y, f) , and X a are the eigenvalues 
of N. 

The equations to be solved are now (75), (76), (82), (90). 

Consider first the modified Vlasov equation (82). The characteristics of 
this equation are defined as the solutions Vj(s,t,v) of the system 

^ = -l^VW^rVsy^VkVtC^ (91) 

with Vj(t,t,v) = vj. Then the solution of (82) is given by f(t,Vi) = 
f (Vj(0,t,Vi)) , and is smooth by standard theory. 
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Now define an iteration as follows: Let /°, a°-, (b )^ be the initial data as 
above, and let = = 0. If smooth iterates a n , b n , m n , w n , f n are 
given for some n, let y n = Lu n , and determine V n+l by solving (91) 
with b n in the rhs, and let f n+1 (t,v) = f°(V n+1 (0,t,v)) which solves the 
n-th Vlasov equation. Now substitute a n , b n , y n , f n+1 into the rhs of (75), 
(76), (90). By Theorem 1 of (Rendall and Schmidt 1991) these linear ODE's 
have a unique smooth solution a n+1 , b n+1 , y n+l on some small time inter- 
val [0, T' n+1 ). Let [0, T n+ \) be the largest time interval on which (det a 11 ) re- 
mains strictly positive. Then by standard theory for regular, linear ODE's, 
the (n+l)th solution can be extended uniquely and smoothly onto [0, T n+ \). 

Let \a\ be the maximum modulus of any component of (a^) , with a simi- 
lar definition for other quantities, and suppose that Vn < N the following 
bounds hold: 

\a n - a°\ < Ai \b n - b°\ < A 2 
\m n \ < A 3 \w n \ < A 4 
(det V 1 )- 1 < A 6 (A 6 > (det b )' 1 ) (92) 
Suppose also that f n (t,v) ^ {v T v) 1 / 2 < A 5 

Now the speed of propagation for the nth Vlasov equation, at the point Vi will 
be bounded by Ci^v) 1 / 2 , with C\ = Ci(A 2 ,Aq). Hence if 
R = diam supp f N (t) , then 

— <dR R< P exp(Cit) (93) 

where Pq = diam supp f°. 
Now (75), (76) give 

\a N - a°\ < C 2 t \b N - b°\ < C 3 t (94) 

and then one gets 

(det fe^r 1 < A % (95) 

for t in a small enough interval [0, T], where T = T(Ai). Since /° is bounded 
one also gets: \H n \ < C^Ai) for n < N. 

From (90) there follows 

vZ(jk)=t-*- fs x «H%-\s)ds (96) 
J o 
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and hence \y& (t)\ < C±t . Also, from the bound on the speed of propa- 
gation, if /° is supported outside r = B\ , then f N is supported outside 
r = Biexp(-Cit). 

Now choose > Po 5 and choose T = T(Ai,Cj) small enough so that 
R < A 5 , and for t G [0, T] f N {t) is supported on {v : <p{v) = 1}, and the 
following bounds hold: 

\a N - a°\ < A 1 \b N - 6° I < A 2 

\m N \ < A 3 \w N \ < A 4 (97) 

Then by induction all iterates exist, are smooth, and are uniformly bounded 
on [0, T]. Also, Vn, f n is supported on : <f>(v) = 1}, and diam supp(/ n ) < 
A 5 . ' 

For the difference between successive iterates one has the following estimates: 

\a n+l - a n \{t) <C f \y n - y 11 " 1 ^) ds (98) 
Jo 

| 6 n+i -b n \(t)<C f \b n - b^Ks) + |y n - y n_1 |(s) (99) 
jo 

_ ^((f) < cj*{(\a n - a 71 ' 1 ] + \b n - b n ^\ + \y n - y n ~ l \){s) 

+ {f (|a"-a- 1 | + |6--6- 1 | + |^-y"- 1 | + ||r +1 -r|| 00 )( S p)#)}d S 

(100) 

Equation (98) comes from (75), (99) comes from (76), and (100) comes from 
(96). 



From the characteristic equation (91) one gets an inequality similar to 
that in (93) and hence the V n are uniformly bounded, above and below 
for « in a compact set with <j){v) = 1 , and t < T, < s < t. Now 
suppose v € {v : <p(v) = 1} PI {v : (v T v) 1 / 2 < ^5}. Then from (91) there 
follows the estimate: 

— -— (s,t,v) < C{\V n+1 - V n \ + \b n - ft"" 1 !} (101) 

CLS CLS 

for t < T, 0<s<t. 
Define the quantity a n (t) by 

a n (t) = sup{|y n+1 - V n \(s,t,v) : < s < t, r < A 5 , 4>{v) = 1} 
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+ \a n+1 - a n \(t) + \b n+1 - b n \(t) + \y n+1 - y n \(t) (102) 

Then 

||/ n+1 -/loo(*)< \\f°\\c^ n (t)) (103) 
Combining the above inequalities gives 

a n {t) < C J* |a n (s) + a n -\s) + j\ a ( s P) + ^(sp)) dp] ds (104) 
Now Gronwall's inequality implies 

a n (t)< C J*{a n - 1 (s) + ^j\ n (sp) + a n ^(sp) dp)} ds (105) 
for t € [0, T], and since the rhs is increasing , this implies 

sup (0 < s < t) a™(s) < C J^{a n -\s) + [ j\a n (sp) + a n ~ l (sp)) dp)}ds 

< Cit(su P(0 < s < t) Q n (s)) + C 2 f sup (0 < p < s) (Q n - 1 (p)) ds (106) 

J 

Put /? n (i) = sup (0 < s < t) a n (s). Then (106) is just 

n (t){l - Cit) < C 2 /"* /^(s) ds (107) 
■/ o 

Choose Ti < T so that (1 - Cit) > \ for i £ [0,Ti]. It follows now that 

(3 n (t) < 2C 2 f I3 n ~ l {s) ds (108) 
J o 

and hence 

P n {t)< C n - 2 \\p 2 \\ 00 t n - 2 /(n-2)l (109) 

Therefore y n , b n , a n , y n are uniformly Cauchy on [0, T±]. It follows 
that a n , b n , y n have continuous limits a(t), b(t), y(t) on [0,Ti], and 
V n (s,t,v) has a limit V(s,i, t>) continuous in s,t. From (75), (76), one gets 
that a, b are C 1 , and from (91) V(s, t, v) is C 1 in s, i, so that a, 6, V in fact 
satisfy these equations. It is then standard (Taylor 1996) that V is C 1 in v. 
If we put f(t,v) = f°(V(0,t,v)), then / is supported on {v : 4>{v) = 1} and 
hence satisfies the Vlasov equation (74) . 

Note that we can differentiate (91), and so V is C 2 along with /. 
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Clearly from (90), (75), (76), (91) we get uniform convergence of the 
time derivatives of all iterates on [e, Ti] for arbitrary < e < T\ , and in- 
ductively it follows that all the sequential limits are smooth on (0, Ti]. 

If we let u = (m,w) = L~ l y as before, then u is continuous on [0,Ti] , 
smooth on (0,Ti] , and satisfies (83) there. Putting k^ = kfj + rriij, Zij = 
Zfj + Wij gives that kij, are continuous on [0,Ti] , smooth on (0, Ti] , 
and satisfy (62), (63) there. These equations can be written 

dr 1 1 

— + -N(a, b, f)x = Gi (a, b, x) + -G 2 (a, b, f) (110) 

where x stands for k^, Z^ , G\ is quadratic in x, G2 accounts for the 
terms 647rG^ p . . . , and Z G (0, T ± ]. 

Now regard a, b,f as known C 1 functions on [0, Ti]. Then (110) implies 

dx 

Z— +Nx = (N- N{Z))x + ZF^Z, x) + F 2 {Z) (111) 
ciZ 

where TV is as in (83), Ti is C 1 in Z and quadratic in x, and F 2 is C 1 . We 
can now differentiate (111) on (0,Ti] to get 

Z ^ + ( I + N )l = ZR ( Z )<1 + s ( z ) ( 112 ) 
where q = and R, S are continuous functions on [0, Ti]. 

It follows by a lemma of (Rendall and Schmidt 1991) that (112) has a 
unique C 1 solution on (0, Ti], and that this extends to a continuous solution 
of the associated integral equation on [0, Ti]. Hence x is C 1 , and (a, b, k, z, f) 
is a classical solution of the conformal Einstein- Vlasov equations. Induc- 
tively one gets that this solution is in fact C°°. 

If two solutions with the same initial data are given, then it is possible 
to derive an inequality for their difference similar to that obtained for the 
quantity (3 n . It follows that the solution constructed above is unique. 

We have therefore proved the following: 

Theorem 6.1. Let G be a 3-dimensional Lie group of some Bianchi type n, 
and let (e l ) a be a basis of left-invariant one-form fields on G . Suppose 
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we are given a smooth function F° on the cotangent bundle of G such 
that F°(x,p) = f°(pi) where pi are the components of p in the frame (6^ a . 
Suppose also that f° is compactly supported, supported outside a neigh- 
bourhood of the origin, and the integral constraint (66) holds. Then there 
exists exactly one smooth Bianchi type n solution {f(t,Pi),(jij(t)) of the 
massless Einstein- Vlasov equations on G x (0, T] with an isotropic singular- 
ity, satisfying f(t,pi) -> f°(pi) as t -> . 

Note that since the initial 3-metric a?- is determined by f° only through the 
integral relation (67) there will be many (anisotropic) /° which give rise to 
the same a^-. But then is determined by the third and fourth moments 
of /°, and so we expect in general to be able to get different values of k®j , 
and thus different 4-geometries, from the same starting metric a^-. This is 
in contrast to the perfect fluid cosmologies of ATI, where it was shown that 
the 4-geometry is in 1-1 correspondence with the 3-geometry on S. 

6.6 FRW data 

The Bianchi types I, V and IX contain, respectively, the k = 0, —1 and +1 
FRW models. We now discuss the question of which f° constitute FRW 
data for the conformal EV equations in these Bianchi types. 

(a) Type I 

The situation here is very simple, as all the C l - k are zero, and one solution of 
the field equations is f(Z,pi) = f°(j)i), g%j = Z 2 ^- , for any admissible f°. 
Hence Theorem 6.1 gives that all Bianchi I solutions of massless EV with 
an isotropic singularity are in fact FRW and the Weyl tensor is always zero. 

(b) Type V 

Here the structure constants are 

C) k = 5/^-5^ (113) 

for some vector a%. 

All type V 3-metrics have constant (negative) curvature, and for f° to 
be FRW data one must certainly have also k^- = 0, since trfc = 0. An 
isotropic /° satisfies this condition by (73), and gives the FRW evolution 
f = f°, (jij = R 2 (Z)a®j, for some function R. It is not clear whether f° must 
be isotropic to give an FRW solution (see (117) below). 
Leaving this question aside for the moment, we consider a related question, 
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namely whether one can pick an anisotropic f° so that the Weyl tensor 
vanishes at E , but the evolution given by Theorem 6.1 is not FRW. We 
calculate that the electric and magnetic parts of the Weyl tensor defined 
relative to V a Z vanish at Z = if and only if 



n 1.0 lm _ n 
u m K l( je^ — u 



kl(k°y j + P ij + \a ij = 



(114) 
(115) 



where is the volume form, A just takes care of the trace and 
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(^0jV2 



flS 



df° (v T k°v)(b )i r v q v s ViV 



3 ^ 



dv t 



<fv+- 



24vr 



_647T_ k 

(a ) 1 / 2 np 



+2C s tr {V 



(v T b°v) 5 / 2 (a°)V2 
-2(k°r(k°y nXijef 

-(k°) ln 

. ln r df v sVq (b°r 



f {v T k v) 2 VtV:j 

(v T Wvfl 2 V 



f° f fxP ^ .., n 2v l v j v k v l v m (b°) m P 

(^TfrO^j [i 6 i v i v kVi + perm ijkl} - - 



{v T b°v) 



~{$ P i v j v kVi + perm ijkl} d v 



+ 



dv t (v T b°v) 3 / 2 
(b°) ln J ^^v^vjVkVi + perm ijkl} d 3 v 



[ f° 

W°) ln J ( v T b o v y v i v 3 v kViv m 



2{b°) m P(v T k°v) 
(v T b°v) 



(k ) mp ) d 3 v 



(116) 



where a =deta o -. 



Now suppose that one has an FRW Einstein- Vlasov metric: a^- = R 2 (Z)a ( -j. 
We then calculate from the field equations that the following condition must 
hold: 



/ 



d n f v iVj :1 



dZ n (v T b°vy/ 2 



d 6 v 



0. 



(117) 



for all n > 1 (it is not clear whether (117) this is sufficient to force isotropy 
of/)). 

From the Vlasov equation (57) we get for n = 1, 2 that conditions (117) 
are respectively: 

Xljk a k = (118) 
WXijkia k a l - 2a iaj - (a fe a fc )a°- = (119) 
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Suppose now that we pick an f° so that (118) and (119) hold, but (117) 
fails to hold for some n > 2 (note that (118) and (119), together with (66) 
and (67), amount to a set of conditions on the first four moments of f° 
only, so that there are many non-isotropic f° satisfying them, while (117) 
for larger and larger n imposes conditions on ever larger moments of /°). 
A calculation then shows that (114) and (115) are satisfied with = 0, so 
that the Weyl tensor vanishes at S. Thus it is possible, with such an f°, to 
generate a non-FRW Einstein- Vlasov cosmology with an initially vanishing 
Weyl tensor. 

(c) Type IX 

Here there exists an invariant basis for which Cj fe = €ijk, and a 3-metric 
hij has constant curvature iff oc 5^ in this basis. So if we let be such 
a metric then kfj = by (73) and the cosmology given by Theorem 6.1 is 
then (k=+l) FRW with f{Z, Pi ) = f°( Pi ), ~ gij = R 2 (Z)a^, for some R(Z) 
regardless of the choice of f°. Thus FRW data in type IX is any admissi- 
ble /° which gives rise to a constant curvature 3-metric. 
As in type V one can ask whether an /° can be chosen so that the Weyl 
tensor vanishes at S, while the 4-geometry departs from FRW shape. Here 
the condition C a b c d(0) = amounts to the following 

k% = (120) 

(95 l e 5/- Xi ; f )s ef = lQ tJ (121) 
where Sij is the trace-free Ricci tensor of a?- and 

Qij = -l^) ln {b°y S Cl s C k n{l {5 k P Xrtl\ 3 ) + ^Xrtk\ 3 ) + S^Xrtkl} (122) 

We have not been able to find a solution of these equations for which Sij 7^ 
and it seems unlikely that there is one. If this is so, then any / which 
makes the initial Weyl tensor vanish is FRW data. 
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